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Abstract 

We elaborate on the fact that quarkonium in hot QCD should not be thought of as a stationary 
bound state in a temperature-dependent real potential, but as a short-lived transient, with 
an exponentially decaying wave function. The reason is the existence of an imaginary part 
in the pertinent static potential, signalling the "disappearance" , due to inelastic scatterings 
with hard particles in the plasma, of the off-shell gluons that bind the quarks together. By 
solving the corresponding Schrodinger equation, we estimate numerically the near-threshold 
spectral functions in scalar, pseudoscalar, vector and axial vector channels, as a function 
of the temperature and of the heavy quark mass. In particular, we point out a subtlety 
in the determination of the scalar channel spectral function and, resolving it to the best of 
our understanding, suggest that at least in the bottomonium case, a resonance peak can be 
observed also in the scalar channel, even though it is strongly suppressed with respect to the 
peak in the vector channel. Finally, we plot the physical dilepton production rate, stressing 
that despite the eventual disappearance of the resonance peak from the corresponding spectral 
function, the quarkonium contribution to the dilepton rate becomes more pronounced with 
increasing temperature, because of the yield from free heavy quarks. 
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1. Introduction 



Assuming the existence of a thermalized medium, with a temperature T, and of heavy quarks, 
with a mass M ^ T, there is a finite probabihty, given by the Boltzmann factor exp(— 2M/r), 
that an on-shell quark and antiquark are generated through thermal fluctuations. They could 
then annihilate, creating an off-shell photon, which may escape from the thermal system, and 
subsequently decay into a dilepton pair (for instance, e~e~^ or fi~fi^). The characteristics 
of the energy distribution of these pairs offer an indirect probe on the strongly interacting 
dynamics taking place within the thermalized system. As a concrete application, properties 
of lepton pairs can be observed in heavy ion collision experiments (see, e.g., refs. [l]), and 
may serve as an indication of whether a thermalized state with a temperature above the 
deconfinement transition was momentarily reached during the evolution [2]. 

Given that various properties of the quarkonium system can be understood in great detail 
at zero temperature [3j, it could be assumed that describing quantitatively the heavy quark- 
antiquark system in a thermalized medium is a relatively simple task. After all, QCD is 
asymptotically free, so the effective coupling should decrease with the temperature, and 
ultimately confinement is lost as well. Somewhat surprisingly, this expectation appears to be 
overly optimistic. In fact, all standard approximation methods develop further systematic 
uncertainties at T > 0. For instance, direct lattice QCD reconstructions of the quarkonium 
spectral function [H El E], which is a quantity determining the dilepton production rate, 
develop the new problem that an analytic continuation is needed from data collected on a 
short Euclidean time interval, to the observable defined in Minkowskian spacetime. Another 
popular class of approaches, so-called potential models [3 [8], suffers from the proliferation 
of many independent non-perturbative definitions of a "static potential" which could be 
measured on the lattice OdD] and inserted into a Schrodinger equation. A recently introduced 
method, the determination of the corresponding observable in strongly coupled = 4 Super- 
Yang- Mills theory [11], also contains unknown systematic errors from the point of view of 
QCD, which cannot be reduced by increasing the temperature, because the QCD coupling 
soon becomes weak |12] . 

The method employed in this paper is resummed weak-coupling perturbation theory. It 
again suffers from novel difficulties at finite temperatures: curing infrared divergences ne- 
cessitates carrying out complicated resummations [13\ 114]. and even though a weak-coupling 
expansion in the QCD coupling constant g can subsequently be defined, it has a strange struc- 
ture, with relative corrections suppressed only by odd powers of g (IS HE], by logarithms like 
g"'ln{l/g) [TTl [18], or by powers of g multiplied by non-perturbative coefficients [T9]-[21j. 
Moreover, even if a number of coefficients were known, the convergence of the series could be 
slow [22] (see, however, refs. [23| [12]). 

Given all these problems, a suitable practical approach at the present date might be to 
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compute the quarkonium spectral function and the dilepton production rate with many dif- 
ferent methods, possessing complementary systematic errors, and to look for a consistent 
pattern, which could then also represent the situation in QCD. It is in this spirit that the 
purpose of the present paper is to pursue the side of resummed perturbative computations. 

The resummed perturbative approach to heavy quarkonium in hot QCD was initiated in 
refs. [2l]-[26], of which the present paper is a direct continuation. In particular, we expand 
and improve on the analysis of ref. ^5]. We consider, first of all, the same observable as in 
ref. [25j (quarkonium spectral function in the vector channel), but discuss more extensively 
the dependence of the result on the temperature and on the heavy quark mass. Second, we 
carry out a new analysis for the spectral function in the scalar channel. This turns out to 
require more advanced numerical techniques than those used in ref. [25|. We also relate the 
pseudoscalar and axial vector spectral functions to the vector and scalar spectral functions. 
Finally, we elaborate on the physics implications of the results in more detail than before, 
both conceptually, i.e. with regard to the picture they suggest for the quarkonium system 
in a deconfined environment, and from the practical point of view, i.e. with regard to the 
dilepton production rate. 



2. General framework 

We start by specifying somewhat more quantitatively the main ideas and equations of the 
resummed perturbative approach. A detailed derivation follows in Sees. O HI while a reader 
only interested in the numerical results could skip directly to Sec. [5] after the present section. 

Let up he a. generic heavy quark field operator in the Heisenberg picture. The basic corre- 
lation function we consider is of the form 

// / 
d3x(|(t,x + 07^W^V'(t,x-0 ^0,-^)j^W'i,(0,+^)) , (2.1) 

where W, W are Wilson lines connecting the adjacent operators, inserted in order to keep the 

Green's function gauge-invariant; the metric is rj^y = diag(H ); and the expectation value 

refers to (...) = 2^~^Tr [exp(— i?/r)(...)], where Z is the partition function, H is the QCD 
Hamiltonian, and T is the temperature. The superscript in CY, refers to the vector channel; 
the subscript refers to the time-ordering in Eq. (j2.ip . We also consider scalar, pseudoscalar 
and axial vector correlators below; their precise definitions are given in Sec. [31 H 

The significance of the Green's function in Eq. (j2.ip is that if we take the limit r, r' 0, 
and subsequently Fourier transform with respect to the time t, then we obtain a function 

^In ref. [24], we set r' — 0, and denoted the correlator by C^{t,v) = C>(t;r, 0). However, with certain 
channels, it will be advantageous to keep r' 7^ 0, because then the singularities from the static potential at 
r = 0, and from the initial condition of the Schrodinger equation at r = r', do not overlap. 
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which is triviaUy related to the heavy quarkonium spectral function, {uj), in the vector 
channel: 

/,^H = i(l-e-^) y""dte^-*C>(t;0,0). (2.2) 

This quantity is physically important, given that the production rate of pairs (with 

a vanishing total spatial momentum = q^- + and a non- vanishing total energy ui = 
E^- + E^+) from a system at a temperature T, is directly proportional to {oj) [27] : 

diV„-„+ 2c^e^ ( 2ml\ / Aml\^ r i/ 1 



d^xd^Q 3(27r)5u;2 \ ' oj"^ J \ ^2 

where we assumed uj > 2m^; e is the electromagnetic coupling; c G (| > ~3 ) is the electric 
charge of the heavy quark; and n-B^uj) = l/[exp(a'/T) — 1] is the Bose distribution function. 

Now, a systematic perturbative determination of the Green's function in Eq. (|2.1|) . and 
of the corresponding spectral function in Eq. (j2.2|) . is quite difficult for energies lv close to 
the quark-antiquark threshold, to ~ 2M. The reason is that in this regime infinitely many 
graphs, particularly so-called ladders, contribute at the same order. A further problem is 
that at finite temperatures, the rungs of the ladders, containing gluons, need to be dressed 
by thermal corrections. 

A way to resum these infinitely many dressed contributions is not to compute the corre- 
lator of Eq. (j2.ip directly, but rather to find a partial differential equation satisfied by this 
correlator, and then to solve this equation numerically. The partial differential equation in 
question is just the Schrodinger equation. Indeed, it is for the sake of being able to write a 
Schrodinger equation that we have introduced r, r' 7^ in Eq. ()2.ip . To be more precise, let 
us consider C> in the limit that the heavy quark mass M is very large. Then, as we will see, 
obeys 

- [2M + y>(t,r) - ^ + (^{-^)] }^>(*'^'^') = ' (2-4) 

with the initial condition 

C^(0;r,r') = -67V,5(3)(r_r') + o(^l:) , (2.5) 

where A'c = 3. The terms specified explicitly in Eqs. (j2.4p . (12. 5p result from a tree-level 
computation; in contrast, the potential denoted by V^{t,r) originates only at 1-loop order. It 
can be defined as the coefficient scaling as 0{M^), after acting on C^(t;r, r') with the time 
derivative idt- The potential V>{t, r) depends, in general, on the temperature; we assume that 
T is parametrically low compared with the heavy quark mass, T ~ {g'^...g)M (cf. Sec. 13. ip . 

Now, it can be argued that in order to be parametrically consistent, the static potential in 
Eq. ()2.4p should be evaluated in the limit t ^ r (cf. Sec. 13. ip . Then it obtains a simple form: 
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in dimensional regularization (cf. Eqs. (4.3), (4.4) of ref . [24j ) . 



limV>(t,r) = — mo H cpimDr) + 0{g*) 

t^oo 47r r 47r 



(2.6) 



where Cp = {N^ — l)/2Nc] is the Debye mass parameter; and the function 




(2.7) 



is finite and strictly increasing, with the limiting values (f)(0) = 0, (j){oo) = 1. 

The first term in Eq. (12. 6p corresponds to twice a thermal mass correction for the heavy 
quarks (cf. the first term inside the square brackets in Eq. (j2.4p ). The second term is a stan- 
dard r-dependent Debye-screened potential. The third term represents an imaginary part: 
its physics is that almost static (off-shell) gluons may disappear due to inelastic scatterings 
with hard particles in the plasma. This is the phenomenon of Landau-damping, well-known 
in plasma physics. As a consequence of the imaginary part, the solution of the Schrodinger 
equation does not lead to a stationary wave function: rather, the bound state decays expo- 
nentially with time, representing a short-lived transient. 

In the following two sections, we discuss the origin of the formulae presented here, and 
their practical evaluation, in some more detail. We also extend the discussion to the other 
channels. We return to the numerical solution of the Schrodinger equation in Sec. [5j 

3. Schrodinger equation and initial conditions 

Our strategy for the derivation of the Schrodinger equation satisfied by the two-point cor- 
relation functions in various channels will be quite straightforward and "modest" her^: we 
first compute the correlation functions in tree-level perturbation theory, and then expand 
in inverse powers of the heavy-quark mass. At this point we can identify the Schrodinger- 
equation and the initial condition for its solution. Subsequently, radiative corrections are 
expected to multiplicatively correct the terms that already appear at tree-level, and to add 
other terms which are allowed by symmetries, even if they would not appear at tree-level; 
the most important of these is the static potential. As long as there is a hierarchy between 
the different physical scales relevant for the problem (cf. Sec. 13. ip . and we are only sensitive 
to perturbative scales, general principles suggest that the system should remain local in the 
presence of radiative corrections, and that a truncation to a certain order is possible. 

3.1. Power counting 

Let us consider the parametric orders of magnitude of the various terms in Eq. (j2.4p , given the 
potential in Eq. (|2.6p . We recall, first of all, that the term 2M plays no role, since it can always 

^ A more systematic approach might follow by generalizing the framework of PNRQCD [28] to finite T. 
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be eliminated through a trivial phase factor (cf. Eq. (j4.ip ). Around the quarkonium peak, 
the time derivative (or energy) is then of the order of the kinetic terms, i.e. dt ~ 9^ /M. If we, 
furthermore, equate kinetic energy with the Coulomb potential energy (assuming mY)r<l, 
cf. below), we are lead to 

a^, ^ i ^ g^M , 9f ~ - ~ . (3.1) 
r t 

An essential question is now to decide how the temperature, T, is to be compared with these 
scales. Let us assume, first of all, that 

T ~ g^M (case 1) . (3.2) 

Then m-£,r ~ gTr ~ <C 1, and Debye screening plays essentially no role yet: we may assume 
the bound state to exist. In this limit, 

Rey> ~ — ~ > Imy> ~ g^TiniDrf ~ /M , (3.3) 

r 

and the imaginary part can indeed be neglected. 

On the other hand, let us increase the temperature to 

T gM (case 2) . (3.4) 

Then Debye screening plays an essential role, m-£,r ~ gTr ~ 1, and we may assume that the 
bound state has melted: indeed, in this limit, 

Re V> ~ — ~ < Imy> ~ g'^T ~ g^M , (3.5) 

r 

so that the imaginary part of the potential, or the width of the state, dominates over the real 
part of the potential, or the binding energy. 

To summarise, the interesting temperature range is T ~ {g'^...g)M. In principle, para- 
metrically consistent analyses in the two limiting cases may require different methods. In 
practice, we would like to have phenomenological access to the whole range; therefore, in 
the present paper we work (implicitly) in the situation where Rey> ~ Imy>, setting us 
somewhere in the middle of the range. For further reference, let us point out that in this 
situation, rV A ~ rmj)A < A, where A is some gauge field component: the variation of the 
infrared gauge fields is parametrically small on the length scales set by the bound state radius. 

3.2. Vector channel 

Denoting 

y'^(x;r)^V^(t,x+|)7^VFv(t,x-0 , (3.6) 
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where x = (t,x), the vector channel correlator we consider is in general of the type 

C^{x; r, r') = (y'^(x; r)V^{0; -r')) • (3.7) 



For simplicity, we have left out hats from the fields in Eq. (j3.6p . as is appropriate once we go 
over to the path integral formulation in Euclidean spacetime. 

Now, even though we will carry out the computation of Eq. (j3.7p within QCD below, 
it will be useful to rewrite the operators considered in the language of NRQCD [29] (for a 
review, see ref. [30]), because this allows to immediately see their scaling with the heavy quark 
mass M, and because this allows to relate various operators to each other in the large-M 
limit. Following ref. [31], we can start by carrying out a Foldy-Wouthuysen transformation. 



where Dj = dj — igAj, Dj = dj + igAj, and we assume a summation over spatial indices, 
j = 1, 2, 3. Afterwards, we go over to the non-relativistic two-component notation by writing 

= ^ ^ , i^^iOK -0^) , (3.9) 
where we already assumed a representation for the Dirac matrices with 

Here fj^ are the Pauli matrices. Furthermore, it is useful to note that in NRQCD, the actions 
for (j) and 9 are of first order in time derivatives; consequently, one of the degrees of freedom 
propagates strictly forward in time, the other strictly backward in time, and a non-zero 
mesonic correlator at t ^ is only obtained from structures like {(p^{...)9 6\.. .)(!)). 

We now find that for , the leading term with the desired structure is 0{1/M) (this term 
is also a total derivative in the limit r ^ 0). Therefore, the correlator C> is dominated by 
the spatial components . At 0{M^), these become 

y*^(x;r) = et(^t^x + 0afcM^(/<(t,x-0 + (^^^ ^ + (i, x - . (3.11) 

To the extent that interactions between the quark and antiquark are spin-independent (this 
is violated only at 0(1/M)), the Pauli-matrices play a trivial role in the two-point correlator 
made out of these operators, yielding eventually Tr [crfcfj;], if and are being correlated. 

We nowproceed to compute the correlator in Eq. (13. 7p at tree-level. We start in Euclidean 
spacetimqfl and after a spatial Fourier transform (for the moment we keep, for generality, the 



^Since both Euclidean and Minkowskian objects appear in this paper, we try to distinguish between them 
by denoting the former with a tilde. In particular, P — (pQf , p) denotes fermionic Euclidean four-momenta, 
while 7^ stand for Euclidean Dirac matrices, satisfying {'y^,^i^} ~ 2S^i,. Any further unspecified conventions 
can be found in ref. 1241. 
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spatial momentum non-zero, q 7^ 0, unlike in Eq. (|2.ip ). whereby 
Cl(r,q;r,r') 



d^xe 



-«qx 



T,x+- T,x- - V' 0,-- W^V 0,+ 



(3.12) 



-if + M -i$ + M 



P2 + M2 52 + M2 



8iV,^ (2vr)3<5(3)(s - p - q)e^(Pof-%)-+K2p+q)-^ pof^'of + P • s + 2M 

7^. (P2 + M2)(52 + ^ 



POf'% 



(P2 + M2)(52 + ^2) 
i!P_,^(P0f-%)-+^(2p+q)-^ 2^0f^'0f + 4 + 4+q + - 



(27r)=^ 



(3.13) 



where pof = 2iTT{n + ^ ) — i/x, n G Z, denotes fermionic Matsubara frequencies (/i is the quark 
chemical potential), and we have introduced the notation 



Ep = VM2 + p2 



The Matsubara sums can be carried out, by making use of 



pI + e^ 



1 

2E 
1 



valid for < r < /3. This yields 



2 L' 



(3.14) 

(3.15) 
(3.16) 



2 X 



(2^)3 

X S^UFiEp + ^JiW{Ep+^ - ^)e('3-)(^p+^p+.) 
+nF{Ep + /i)nF(i?p+q + ^)e(/3--)^p+-^p+.+/5A' 



-q2 + 2M2 + {Ep + Ep+^f 

EpEp+q 
q2 _ 2M2 - (£;p - ^p+q)2 



q2 _ 2M2 - (^p - Ep+^ 
EpEpj^c^ 



i2n 



+ 



+ 



+ 



-q2 + 2M2 + {Ep + Pp+q) 
£^p^p+q 



(3.17) 



In order to simplify the expression somewhat, we note that once we go over into the spectral 
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functior0, and restrict to frequencies (energies) around the quark-antiquark threshold, \uj — 
2M\ <C M, then only the first of the structures in Eq. (j3.17p contributes. Second, close 
enough to the threshold, the J- function expressing energy-conservation, 5{ijO — Ep — E'p+q), 
forces the loop momentum p to be small, |p| <^ M. We also assume the external momentum 
to be small, |q| <^ M. Under these circumstances, we can expand 



and the relevant part of C^{t, q; r, r') becomes 
~ -6iVc 



M + 



|p_+q|_ 

2M 



(3.18) 



^ ^(2p+q).i:^-.[2Af+ ^P^+^pW +0(^)] 



l + O 



1 



(3.19) 



Here we have also omitted effects of relative order exp(— [M it /i]/T), by keeping only the 
leading terms in the exponentials. We note that after these simplifications, all dependence 
on the temperature and on the chemical potential has disappeared from the tree-level result. 
The real-time object we are ultimately interested in, is the analytic continuation 



Noting from Eq. ([3l9]) that 



C7^(t,q;r,r')=C|(it,q;r,r'). 



iVr P + 



(3.20) 



(3.21) 



the dependence on r and t in the exponential amounts to satisfying the Schrodinger equation 

^2 



idt 



2M + 



AM 



M 



+ }c^(t,q;r,r') = 0. (3.22) 



The initial condition for the solution is obtained by setting t = in Eq. (j3.19p (after use of 
Eq. dS^Ql)): we find 



C^(0,q;r,r') 



-6N, 



d3 



(27r)3 



P_gi(2p+q)- 



r — r 

— +0 



1 
M2 



r') + 



1 

w 



(3.23) 



Eqs. ([3:22]) . ([3:23|) justify Eqs. ([23]), ([23]) for the vector channel in the free limit. 

For future reference, let us also compute (uj) explicitly (general expressions for free 
spectral functions can be found in refs. [38]). Eq. (j3.19|) (after r it) already shows the 



*Take first a Fourier transform, Csii^iy) ~ /p''dr e'"b'^C'B(r), wfiere LJ■^^ is a bosonic Matsubara frequency; 
then carry out tlie analytic continuation p{u;) — ^ [C'b(— + iO^]) — C'b(— — jO+])]. A typical term in 
Ce{t), of the form exp(Air + A2(/3 - r)), becomes p{lo) = -7r(e'"^i - e'^^^)S{Lu + Ai - A2). 
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solution of Eqs. (j3.22p . (j3.23p . and we can then directly remove the point-splitting, setting 
r, r' = 0. Shifting p ^ p — q/2; taking the steps in footnote and ignoring exponentially 
small terms and terms suppressed by 0(1/M^), we find 

- I ^5{.' - g) = -^^(.')mI (.0^ , (3.24) 



where 



uj' = to 



(3.25) 



In the following, we will often for simplicity restrict to q = (like already in Eq. (j2.4p ). but 
we can now observe from Eqs. (j3.24p . (j3.25p that the main effect of a non-zero q 7^ is simply 
to shift the threshold location 2M by the center-of-mass kinetic energy q^/4M. 

The analysis so far has been at tree-level. As argued in refs. |24l I25j . however, the essen- 
tial (temperature and w-dependent) 1-loop corrections can be taken into account simply by 
inserting the potential Vy{oo,r), given in Eq. (j2.6p . into Eq. (j3.22p . There are of course also 
other loop corrections, related for instance to the renormalization and definition of M as a 
pole mass, and the overall normalization of the non-relativistic vector current in Eq. (j3.1ip : 
these corrections are in fact known to high order at zero temperature [32| [33] Fl but are not 
essential at our current resolution, so we mostly omit them here. 



3.3. Scalar channel 

Denoting 

5(x;r) EEV;(t,x+|)l^^(t,x-|) , (3.26) 
the scalar channel correlator we consider is of the type 

C|(x; r, r') = (^S{x; r)5(0; -r')) . (3.27) 

The correlator C>(x; 0, 0) is not directly physical, but it does have the appropriate quantum 
numbers to give a contribution to the three-particle production rate qq /i~/Lt"'"7, i.e. a 
lepton-antilepton pair together with an on-shell photon. Moreover, it is frequently measured 
on the lattice, which will be our most direct reference point. We will ignore the issue of overall 
(re) normalization in the following, and concentrate on the shape of the spectral function 
(meaning its w-dependence in frequency space, or its t-dependence in coordinate space). 

It is again helpful to write S{x; r) with the NRQCD notation. The steps in Eqs. (|3.8p . (|3.9p 
indicate that at 0{M^), S = 9^6 — (fi^cj), which does not lead to any non-trivial t-dependence. 

^To 1-loop order, M = my^{mj;is)i^ + d^CF/^TT^), V^nQcu{x;0) = V^cd{^;0){1 + g^Cr /2tt^). 

^It may be noted, for instance, that the scalar density requires renormalization, unlike the vector current. 
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The leading non-trivial term reads 



S{x;r) = ... + — 
^ ' 2M 



where = VKD/t, x - r/2) - D/t, x + r/2) VF. 

To simplify Eq. (j3.28p a bit, let us for now assume that the gauge fields are perturbative, 
so that the Wilson line can be approximated by the first term in its expansion; and that 
their variation is slow on the scale set by |r|, as argued in Sec. 13.11 (in any case, |r| is 
taken to be zero at the end). Then we may write ~ 1 + igr ■ A(t,x), Dj{t,x — r/2) ~ 
'dj-igAj{t,x.) +igr ■VAj{t,x)/2,Dj{t,x + r/2) 'dj+ igAj{t,x.) + igr ■ VAj{t,x)/2. We 
now note that 

(t, X + I) ^(t, X - -2^ {^t X + 1^0 X - £) } . (3.29) 

Therefore, to leading order in the large-M expansion, and at least to some order in the 
weak-coupling expansion, we can identify 

5(x;r)~-^Vr-V(x;r) , (3.30) 

where the components of V are given in Eq. (13. lip . 

The relation between the vector and scalar channel correlators can be pushed one step 
further, if we consider directly the correlators, Eqs. (13. 7p and (I3.27p . To leading order in the 
large-M expansion, Eqs. ()3.1ip and (I3.30p imply that 

C|(x;r,r') = /^(x; r)5(0; -r') 



1 ^ 

T2 5](Vr)fc(Vr')/(^fc(x;r)y/(0;-r' 



~ M2 

k,l=l 

3 3 

= (Vr)fc(VrO/ SMj2{Vj{x;r)V,{0;-r': 

k,i=i j=i 

= -3i2^r-Vr'^(y^(x;r)lS(0;-r'; 



^ rVr- Vr' C^{x;r,r') . (3.31) 



3M2 

We will be making use of this important relation later on. 

We now return to full QCD, and outline the computation of the 2-point scalar density cor- 
relator in Eq. (I3.27|) at tree-level, again taking a spatial Fourier transfrom and, for generality, 
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keeping track of a non-zero spatial momentum q 7^ for the moment. Then, 
C|(r, q;r,r') 



d'^xe-*'i-"(^ T,x+- W' r,x-- 0,-- V 0,+ 



(3.32) 



-if +M -i$ +M' 



d3 



(2vr)^ 



(3.33) 



Making use of Eqs. p.l5p . (j3.16p . this can be rewritten as 



C|(r, q;r,r') 



-Nr. 



P*(2p+q)- 



2 X 



X |nF(i^p + ^i)np{Ep+^ - /,)e(/5--)(^p+^p+'>) 
+nF(^p + fi)nF{Ep+q + ^)e(/3-)^P+-£^P+<i+/3M 
+nF(£;p - ^)nF(-E;p+q - ^)e(/3-^)^p+q+^^p~/3M 



i2n 



+ 



q2 + 4M2 - (£;p + Ep+^ 

2-E'pi?p+q 

-q2 - 4m2 + (£;p - £;p+q)2 

2£'p-E'p+q 
-q2 - 4M2 + (^p - ^p+q)2 
2-E'p-Ep+q 



+ 



+ 



+nF(£;p - /i)nF(£;p+q + /i)e^(^p+^p+<i) 



q2 + 4m2 - (^p + Ep+^) 

2-Ep-E'p+q 



(3.34) 



With the same considerations as between Eqs. p.l7p and (j3.19p . the interesting part of Cj 
can be approximated as 



C|(r, q;r,r') ^ ^ /^e^{2p+q)-^-[2M+2E!±||^+o(^)] 



1 
M2 



2M2 J (27r) 
4p2 + 4p • q + q2 + O 



(3.35) 



Note again that after these simpUfications, all dependence on the temperature and on the 
chemical potential has disappeared from the tree-level result. 

The exponential in Eq. (|3.35p is the same as in Eq. (j3.19p . whereby C> obeys the same 
Schrodinger equation as C^, Eq. (|3.22p . The initial condition is different, however: setting 
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t = in Eq. ()3.35p (after t —>■ it), we find 



Cf(0,q;r,r') 



2N, 



M2 

2iVe 

MS 



(27r)3' 
V2 5(3)(r_: 



P*{2p+q) 



r — r 

~2~ 



(3.36) 



Tliis agrees, of course, with what can be deduced from Eqs. (j3.23p . (j3.3ip . We note that all 
dependence on the external momentum q again only appears as a part of the center-of-mass 
energy 2M + q^/AM, inside Eq. ([3:22]) . 

For future reference, let us finally determine the spectral function, p'^[uj). Eq. (j3.35p (after 
r — > it) already shows the solution of Eqs. (|3.22p . (j3.36p . and we can then directly remove 
the point-splitting, setting r, r' = 0. Shifting p ^ p — q/2; taking the steps in footnote [H 
and ignoring exponentially small terms, we find 



p'ico) 



M2 



f <^ P 2.( I P 

7 (^p "M 



— e{J)M'2{j) 

27r 



(3.37) 



where uj' is from Eq. (|3.25p . 

The analysis so far has been at tree-level. As discussed above Eq. (|3.29p . the relation 
in Eq. (j3.3ip is more general, however. Therefore, we can extract a beyond-the-leading order 
by simply applying Eq. (j3.3ip to a beyond-the-leading order p^ . 



3.4. Other channels 



In Sees. 13.21 13.31 we have discussed the correlators in the vector and scalar channels. Let us 
now show that in the limit of a large quark mass, the correlators in the pseudoscalar and 
axial vector channels are to a good approximation equivalent to either of these two. 
We note, first of all, that in the basis of Eq. (j3.10p . the matrix 75 becomes 



75 




(3.38) 



Thereby the pseudoscalar density becomes 
P{x-r) 



V'(t,x + 0i75Ty^(t,x- ^ 



(3.39) 



t,.+ '^we{t,. I 



(3.40) 



where again only structures of the type 6^(j) and (jrO have been kept. The non-trivial two-point 
correlator comes from the cross-term between the two structures in Eq. (j3.40p . and ignoring 
the spin-dependent corrections of ©(l/M), a comparison with Eq. (|3.1ip then shows that 

1 



C>(a;;r,r 



C> (x; r, r 



(3.41) 
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where C^(x;r,r') = (P(x;r)P(0; -r')), and is defined in Eq. ([STTl) . 
The axial vector, on the other hand, can be defined as 



yl^^Cx; r) = ( t, X + I) 757^^ VFV (t, X - I 



(3.42) 



In the case of V^, we found that the dominant contribution is given by the spatial components, 
but for the axial vector, the roles have interchanged: the leading term is 



A\x- r) = - [e^ (t, X + W<p(t, X - ^) + 
Comparing with Eq. (j3.40p . we find 



<,x + |),w(t,x-| 



(x;r,r') = r)A°(0; -r')) ~ C^(x;r,r') ~ --C^(x;r,r') . 



(3.43) 



In lattice studies, however, attention is sometimes restricted to the spatial components A^; 
repeating the previous steps, we find 



+ 



t,x + ^yiwe{t,x-^ 



. (3.44) 



The first term is a total derivative, and the second term has a structure close to that in 
Eq. (j3.1ip . given that only the crossterm contributes in a correlation function. Therefore, 
paralleling the argument in Eq. (|3.3ip . we find 



M2 



dr^dj 



J / d3x(y'(x;r)y''(0;-r')) 



3M2 



^klm^kV' 



^6w y"d3xC^(x;r,r') 
2 j d3xCf(x;r,r') . (3.45) 



To summarize, Eqs. (|3.4ip . (|3.43p and (|3.45p show that the pseudoscalar and axial correlators 
do not lead to any qualitatively new structures. 



4. Method to construct the spectral functions 

In the previous section, we have set up the Schrodinger equation and initial conditions satisfied 
by the vector channel correlator C^, and shown that the corresponding correlators in the 
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other channels can be obtained from through various relations. The aim now is to extract 
the spectral functions corresponding to these correlators. 

To achieve this goal, it is useful to convert the time-dependent Schrodinger equation directly 
to frequency space. Let 

^Pit■r,r')^e'^'''C^{t,q■,r,r'), (4.1) 



and 



X{t; r, r') = e^''^'C^{t, q; r, r') ^^^(t; r, r') 



The corresponding frequency representations are defined by 



(4.2) 



and the spectral functions are then obtained from (cf. Eq. (12.2 



i^iu;'; r, r') ^ / dt e'^'' ^(t; r, r') , x{uj'- r, r') ^ / dt e^-'* x(t; r, r') , (4.3) 



(4.4) 
(4.5) 



lim - V'(w';r,r') , 

r,r'— >0 Z 
r,r'^0 Z 



where uj' is from Eq. (|3.25p and we have omitted exponentially small corrections. 

We now recall from ref. [Mj that the imaginary part of y>(t, r) (Eq. (j2.6p ) is odd in t — > —t. 
Furthermore, we recall from Sec. 13. II that a consistent perturbative solution allows (or, to be 
more precise, demands) considering the limit |t| ^ r. Denoting 



y>(r)= lim y>(t,r), 

the equations to be solved (Eq. (j2.4p ) then read 

H — i\lm.Vy{r)\ 'ip{t;r,r') = idttpit; r,r') 
H + i\lm.Vy{r)\ ip{t;r,r') = idtipit; r,r') 



t>0 , 
t<0 , 



(4.6) 

(4.7) 
(4.8) 



where we indicated explicitly that the imaginary part is negative for t — > +oo [24 ^ I26 j . and 
defined a Hermitean differential operator H through 



H 



(4.9) 



Since the effective Hamiltonian is time-independent both for t < and for t > 0, we can 
formally solve Eqs. <SM- 



V^(t;r,r') 



-iHt-\lmV>(r)\t 
-iHt+\lmV>{r)\t 



V'(0;r,r') , t>0 
V'(0;r,r') , t<0 



(4.10) 
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where, according to Eqs. (|3.23p . (|4.ip . 

V'(0;r,r') = -6Nc6^^\r - r') . 
Taking a Fourier-transform, we get 



(4.11) 



dte^'^'V(t;r,r') 

oo 

iJ -iH - |Imy>(r)| 
iJ - iH +\lmVy (r) | 
uj' — H — i\ Im Vy{r) 



M^'-H)-\lraV>{:r)\t 



Jt[u' -H)+\lmV>{r)\t 



+ 



V^(0;r,r') 



J - H + i\\m.V>{r)\ }>V40;r,r'). 

(4.12) 



To give a concrete meaning to the inverses in Eq. (j4.12p . we define a function ^'(cj'; r, r') as 
the solution of the equation 



r — r 



J - H + i\luiV>{r)\ ^'(a;';r,r') 
Then the result of Eq. (j4.12p can be rewritten as 

^{J; r, r') = -2 Im ^{J; r, r') . 
According to Eqs. (|4.2p . (j4.4p . (j4.5p . the spectral functions are now obtained from 



P (^) 



p^{ijj') ~ lim 



lim Im 

r,r'— ►O 



^{uj'\ r, r') 



1 



o3M2 



Im 



Vr • Vr'^'(a;';r,r') 



(4.13) 

(4.14) 

(4.15) 
(4.16) 



To summarize, we have reduced the determination of the spectral functions to the solution 
of a time-independent inhomogeneous Schrodinger equation, Eq. (I4.13p . 
As the next step, following ref. [36], we introduce the ansatz 



oo I 



1=0 m=-l 



(4.17) 



Here Yjm are the spherical harmonics, normalised as /dfi y5*^(r2)yj'm'(^) = ^U'Smm', where 
dQ = dcos 6 d(/), and satisfying 

J2 Yhn{^')yimi^) = - = HCOS 9 - cos 9') (5((/. - (f)') . (4.18) 

Im 
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The J-function can be written as 

(^(3)(r - r') = -^6{r - r')6{n - n') 
whereby Eq. (|4.13p becomes 

Hr + i\lmVy{r)\ gi{J]r,r')-- 



UJ 



with 



(4.19) 



(4.20) 



(4.21) 



The remaining goal is to reduce the problem to the solution of the homogeneous equation. 
Following refs. [SGlfH], we introduce the ansatz 
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(4.22) 



where g^^ is a solution of the homogeneous equation regular at zero; gi^ is a solution of the 
homogeneous equation regular at infinity; and r< = min(r, r'),r> = max(r, r'). 



Obviously, the function gi is symmetric in r <-> r', and continuous at r 
well-known form of the solution gl^ , it must thus behave as 



r . 



Given the 



5i ~ r - + 0(r'+2)] [(r')'+^ + 0((r')'+')] (4.23) 
at small r, r'. For the vector channel spectral function, Eqs. (j4.15|) . (|4.17|) now imply that 



lim ^ 

r.r'^o 47rrr' 



Im 



go{uj';r, r') 



(4.24) 



i.e. that only the S-wave (/ = 0) solution of the homogeneous part of Eq. (|4.2U|) contributes. 

Consider then the scalar channel. According to Eq. (|4.16|) . the scalar channel spectral 
function can be extracted from the same function ^ as the vector channel one, by taking two 
derivatives and then extrapolating r, r' — > 0. Inspecting Eq. (|4.23p . we see that we at least get 
a contribution from the P-wave {I = 1). However, according to Eq. (j4.23p . it is also possible 
to get a contribution from the subleading S-wave terms, go ^ [r + 0{r'^)][r' + 0{{r')^)]. As 
far as we can see, this contribution was omitted in ref. [8|. 

We relegate a more detailed discussion on how to write the solutions for the spectral 
functions , to Appendix A, given that the further steps are quite technical in nature, 
and give here just the final formulae. Introducing the dimensionless variables g = raM and 
a = g'^Cp /4:7r, the vector channel spectral function from Eq. ()4.24p can be simplified to 



M2 



lim 



dglm 



1 



(4.25) 



16 



while the scalar channel spectral function becomes 



lim 



dg Im 



+ 



36 



mo)? mQ)] 



9l(e)=e2-e3/4+.. 



(4.26) 



We remark that because of the factor 36, the first (S-wave) term is numerically subdominant 
in Eq. (14.26p . and would be totally negligible, were it not for the fact that is does lead to a 
resonance peak, unlike the second term. 



5. Numerical results 

In the previous section, we have reduced the numerical determination of the vector and scalar 
channel spectral functions to Eqs. ()4.25p . (I4.26|) . respectively. In these equations the functions 
5^, / = 0, 1, denote the regular solutions of the homogeneous part of Eq. (I4.20p . 

J - Hr + i\liaV>{r)\\g^^ = {) , (5.1) 

where is from Eq. (j4.2ip . Further details can be found in Appendix A. 

In practice, the procedure of determining , starts from some small value, g = 6, 
with for instance 5 = 10~^, at which point we impose as initial conditions the properties of 
the regular solutions at small g, g%{5) = 5 — 5"^ /2 + ... , g\{5) = 5^ - (5^/4 + ... . We then 
integrate Eq. ()5.ip towards larger g, constructing simultaneously the quantities in Eqs. (j4.25p . 
(j4.26p . After a while, g^^{g) and g\{g) start to grow rapidly and the integrals in Eqs. (j4.25p . 
(j4.26p settle to their asymptotic values. Subsequently, we check that the results obtained are 
independent of the starting point 5. The numerics is straightforward and poses no problems. 

Apart from the pole mass M, the solution depends on what is plugged in for g^ and mD- 
We employ here simple analytic expressions that can be extracted from Ref. |37j . 

^'" 91n(9.082r/A^) ' " 31n(7.547r/A^) ' ^r iV. = iV, = 3, ^ = . (5.2) 



We also fix Aj^ ~ 300 MeV; for the uncertainties related to this, see Fig. 2 of ref. [2^. 

The results for -p^ ju? (Eq. (OSl) divided by -ij?lM'^) are shown in Figs. [H [21 and those 
for p^ jij? (Eq. (|06D divided by uP'jM'^) in Figs. El d The results are given in a range of uj 
where relativistic corrections, i.e. terms of higher order in a Taylor expansion in {lo — 2M) /M , 
are estimated to be at most at the 10% level. We show a scan of mass values, given that the 
inherent theoretical uncertainties of the charm and bottom pole masses are several hundred 
MeV (for a pedagogic discussion, see ref. [39]), and that in lattice simulations there are 
further uncertainties, related to scale setting etc, which make it difficult to sit precisely at 
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M = 2 GeV 



M = 4 GeV 



M = 6 GeV 




Figure 1: The resummed perturbative vector channel spectral function p^(a;), in units of w^, 
in the non-relativistic regime, (w — 2M)/M <^ 1, for M = 2,4,6 GeV (from left to right). 
To the order considered, M is the heavy quark pole mass. Note that for better visibility, the 
axis ranges are different in the rightmost figure. 

the physical point. As far as the other channels are concerned, we recall from Eqs. (j3.4ip . 
(fMjl . (fOHjl that 

p^^-lp^ P^'^-Ip""; p^^2pS. (5.3) 
5.1. Comparison with lattice 

As of today, lattice reconstructions of the spectral functions in various channels [H El [6] 
suffer from significant uncertainties. Apart from the usual problems, it may be mentioned 
that the Compton wavelength associated with the heavy quarks tends to be of the order of 
the lattice spacing, so that we may expect even more significant discretization artifacts than 
in the usual quenched or 2+1 light flavour simulations; and that the analytic continuation 
from Euclidean lattice data to the Minkowskian spectral function necessarily involves model 
input, whose uncertainties are difficult to quantify. Nevertheless, it has been claimed that 
the latter types of uncertainties may be under reasonable control from a practical point of 
view UQ]. The most recent lattice results in this spirit can be found in refs. [5l[6]. 

It has become fashionable recently not to compare directly the spectral functions, but 
the Euclidean correlators for which direct lattice data exists. Though this removes the un- 
certainties related to the analytic continuation, it also comes with a heavy price: most of 
the structure in a Euclidean correlator is determined by values of uj far from the threshold, 
u) <C 2M or w ^ 2M, so that the actual physics we are interested in tends to be hidden in 
tiny effects somewhere in the middle of the Euclidean time interval. For this reason, we do 
not consider Euclidean correlators to be as interesting as the spectral functions, and touch 



18 



T = 250 MeV 



T = 350 MeV 



T = 450 MeV 
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IV, 



0.0 
(B/M - 2.0 



1.4 

1.2 

1.0 
j 0.8 
■-0.6 

0.4 

0.2 

0.0 K-^ 
-0.4 



M = 1 GcV 
M = 2 GcV 
M = 3 GcV 
M = 4 GcV 
M = 5 GcV 
M = 6 GcV 
500 McV « T « M 



0.0 
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Figure 2: The resummed perturbative vector channel spectral function p (lo), in units of 
w^, in the non-relativistic regime, - 2M)/M < 1, for T = 250,350,450 MeV (from left 
to right). To the order considered, M is the heavy quark pole mass. Note that for better 
visibility, the axis ranges are different in the leftmost figure. 

only the latter in the following. 

Most of the lattice data exists for the charmonium case. The temperatures where the 
charmonium peak disappears from the spectral function are rather low, however; in fact they 
are in a regime where our analysis is probably not yet justified. Assuming the charmonium 
pole mass to be in the range M ~ (1.5... 2.0) GeV, we nevertheless observe from Fig. [Deleft) 
that at T ~ 250 MeV a certain "enhancement" can still be seen in the vector (and thus, 
in the pseudoscalar) channel. This then disappears at higher temperatures. In contrast, 
in the scalar channel, Fig. [3l|left), there is practically no structure. These observations are 
certainly not in conflict with the lattice results of refs. [5l[6]. Furthermore, we may note that 
the absolute magnitudes of pv and ps in Figs. [Deleft), [Sl^left) are qualitatively in a similar 
relation to each other as the spectral functions measured on the lattice: the difference of 
about an order of magnitude is due to the 1/M^-suppression in the scalar case. At the same 
time, it needs to be kept in mind that in the scalar case the operators require renormalization, 
and that we have in any case not computed radiative corrections to the absolute magnitudes 
of the spectral functions, so that the comparison cannot be taken too seriously. 

Data for the bottomonium case, where our predictions should be more reliable, can be found 
in ref. [5]. There is again an inherent uncertainty of several hundred MeV in the bottom quark 
pole mass, but realistic values are presumable in the range M ~ (4. 5... 5.0) GeV. According to 
Figs. [H [2] (middle to right), there is now a clear peak in the vector channel spectral function, 
up to a temperature of perhaps 500 MeV. In the scalar channel case. Figs. [3l [4] (middle to 
right), the structure is much less pronounced, but a tiny enhancement can be observed up to 
a temperature of about 400 MeV. These results are qualitatively in better agreement with 
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Figure 3: The resummed perturbative scalar channel spectral function p^{uj), in units of w^, 
in the non-relativistic regime, (w — 2M)/M <^ 1, for M = 2,4,6 GeV (from left to right). 
To the order considered, M is the heavy quark pole mass. Note that for better visibility, the 
range of x-axis is different in the rightmost figure. 

the lattice data in ref. [5j than the potential model results of ref. [8], where no peak was 
found in the scalar channel case; as we have explained in Sec. [U the discrepancy can be 
traced back to a difference in the reconstruction of the spectral function from a Schrodinger 
equation. Nevertheless, in practice, it should again be stressed that systematic uncertainties 
of the lattice data are certainly too large to make a quantitative comparison. 

5.2. Dilepton rate 

Apart from the spectral functions, it is interesting to plot also the physical observable, the 
dilepton production rate given in Eq. (j2.3p . This is shown in Fig. [5j The significant difference 
with respect to the vector channel spectral function is the existence of the Boltzmann factor 
(or, to be more precise, Bose-Einstein factor) in Eq. (j2.3p . Obviously, for a fixed frequency 
around the threshold, uj 2M, the Boltzmann factor exp(— a;/T) exp(— 2M/T) introduces 
a strong dependence of the dilepton rate on the temperature or, for a given temperature, 
on the mass. The exponential boosts the rate at high temperatures, and makes it decrease 
rapidly at large frequencies. Thereby the dilepton rate shows a much stronger resonance-like 
behaviour than the spectral function. Fig. [H In particular, some kind of a peak structure 
remains visible in the dilepton rate in Fig. [5] even at temperatures which are so high that 
there is only a smooth step-like behaviour visible in the spectral function in Fig. [H 
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Figure 4: The resummed perturbative scalar channel spectral function /? (u), in units of 
in the non-relativistic regime, (w - 2M)/M < 1, for T = 250,350,450 MeV (from left 
to right). To the order considered, M is the heavy quark pole mass. Note that for better 
visibility, the range of x-axis is different in the leftmost figure. 

6. Physical picture of heavy quarkonium in a thermal plasma 

Conceptually, the most important difference between our analysis and traditional potential 
models [71|8] is the existence of an imaginary part in the static potential, Eq. p. 60 . Physically, 
the imaginary part implies that quarkonium at high temperatures should not be thought of 
as a stationary state. Rather, the norm of its wave function decays exponentially with 
(Minkowski) time. This is due to the fact that, apart from experiencing Debye screening, 
there is also a finite probability for the off-shell gluons binding the two quarks to disappear, 
due to Landau damping, i.e. inelastic scatterings with hard particles in the plasma. Once 
T ~ gM, the imaginary part is in fact parametrically larger than the binding energy (cf. 
Sec. 13. ip . At the same time, for low enough temperatures, T ~ g'^M, the imaginary part 
plays a subdominant role (cf. Sec. 13. ip . 

It may be useful to remark that if, on the contrary, one goes to a Euclidean lattice, then a 
non-zero wave function can be defined at any finite value of the "imaginary time" coordinate 
T, < T < f3. Introducing also gauge-fixing, such wave functions have been measured 
with Monte Carlo simulations in ref. [33] (for a recent review, see ref. [35]). With regard 
to the discussion above, the physical significance of such wave functions for Minkowski-time 
observables is not obvious; hence we do not discuss them here. 
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Figure 5: The physical dilepton production rate, Eq. ()2.3p . from charmonium (left) and 
bottomonium (right), as a function of the energy, for various temperatures. The mass M 
corresponds to the pole mass, and is subject to uncertainties of several hundred MeV; we use 
the intervals 1.5. ..2.0 GeV and 4. 5. ..5.0 GeV to illustrate the magnitude of the corresponding 
error bands. The low mass corresponds to the upper edge of each error band. 

7. Conclusions 

The purpose of this paper has been to experiment, as generally as possible, with the resummed 
perturbative framework that was introduced in refs. [Ml EH], in order to offer one more handle 
on the properties of heavy quarkonium in hot QCD, thus supplementing the traditional 
approaches based on potential models and on lattice QCD. 

The key ingredient of our approach is a careful definition of a finite-temperature real-time 
static potential that can be inserted into a Schrodinger equation obeyed by certain heavy 
quarkonium Green's functions. The potential in question, denoted by \\mt^ooy>{t^f)-, has 
both a real and an imaginary part (cf. Eq. (|2.6p ). An important conceptual consequence 
from the existence of an imaginary part is that heavy quarkonium should not be thought of 
as a stationary state at high temperatures, but as a short-lived transient, with the quark and 
antiquark binding together only for a brief moment before unattaching again. 

On the more technical level we have noted that, in terms of Eq. ()4.17p . the vector channel 
spectral function gets a contribution only from the S-wave, / = 0, while the scalar channel 
spectral function gets a contribution both from the S-wave and P-wave, I = 0, 1. Here we 
differ from the potential model analysis in ref. [8] where, as far as we can see, only / = 1 
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was considered for the scalar channel. The reason for the difference is discussed at the end 
of Sec. m The difference is significant, since the S-wave contribution introduces a small 
reasonance peak to the scalar channel spectral function as well. 

The phenomenological pattern we find for the spectral functions within this framework 
is not too different from indications from lattice QCD: scalar channel charmonium displays 
practically no reasonance peak above a temperature of 200 MeV; vector channel charmonium 
has some peak-like structure up to a temperature of about 300 MeV; scalar channel bottomo- 
nium is again weakly bound but does show a small enhancement up to a temperature of about 
400 MeV; vector channel bottomonium can support a resonance peak up to a temperature 
of about 500 MeV. (Because of unknown higher order corrections, these numbers are subject 
to uncertainties of several tens of MeV.) 

At the same time, we stress that in the physical dilepton rate. Fig. [5l the quarkonium peak 
always becomes more pronounced with increasing temperature, irrespective of the disappear- 
ance of the resonance structure from the spectral function. This boost is due to an interplay 
of the free quark continuum in the spectral function, and the Boltzmann factor exp(— u/T). 

There are a few directions in which our work could be extended, in order to go beyond 
a purely perturbative approach. In particular, the imaginary part of the real-time static 
potential has been measured with classical lattice gauge theory simulations in ref. [26j, and 
could thus to some extent be used in a non-perturbative setting. Hopefully, the real part 
of our static potential could also be related to quantities that are measurable with lattice 
Monte Carlo methods, thereby allowing us to probe more reliably the phenomenologically 
interesting temperature regime around a few hundred MeV. 
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Appendix A. Numerical method for finding the spectral functions 

In this appendix we provide details concerning the numerical method that we have used for 
determining the vector and scalar channel spectral functions. The basic approach is from 
ref. [36], where it was applied for the vector channel at zero temperature; the method was 
extended to the scalar channel case in ref. [8j. Our presentation is rather close to that in 
ref. [8], but we choose to spell out the details anew due to the fact that, as already mentioned 
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in Sec. HI we find one additional term in the scalar channel case. Furthermore, the existence 
of an imaginary part in our static potential simplifies certain points of the analysis. We 
should point out that the method presented here appears to be numerically superior to that 
introduced for the vector channel in ref. 1251. 



A.l. Vector channel 

We proceed with the evaluation of Eq. (|4.24|) . Given the ansatz in Eq. (j4.22p . it remains to 
determine A,g^^,gi^, and then to extrapolate r, r' — > 0. We thus need to know, in particular, 
the asymptotic behaviours of the functions (7^ , gl^ near the origin. Let g[ and g\ be the 
solutions regular and irregular around the origin, respectively: 

00 

g\ = r'+i5]a„r"«aor'+S (A.l) 

n=0 



al = aiir) / dr' 



1 1 r 



-I 



]2 ao 2/ + 1 



(A.2) 



We may then choose 



9'<{r) = giir), (A.3) 
9>(.r) = g\{r)+B'g[{r) , (A.4) 

where the coefficient is defined such as to guarantee the regularity of gi^ (r) at infinity, 

i?' = - hm = - /""dr' -i-- . (A.5) 

Combining Eqs. ()A.2p . ()A.4p . (|A.5p . we can write 

Let us next compute the coefficient A in Eq. ()4.22p . Integrating both sides of Eq. (|4.20p 
with /^^^J^Q^+dr (...), yields 



6NcM 

g'^{r')dg'-^{r') /dr' — g'-^{r')dg'^{r') / dr' 



Involving a Wronskian, this expression is independent of the position r' at which it is evalu- 
ated, so we can do this at small r' . Then we can use the asymptotic forms from Eqs. (jA.ip . 
(jX2|) . to find that 

A = -GNcM . (A.8) 
Note that this expression is independent of /. 
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We finally take the limit r, r' 0, while keeping r < r', so that r< = r, r> = r'. Inserting 
Eqs. g^lD, (!Al]) . and (IXSl) into Eq. g^D, yields 



6iVcM , 1 ^ 
— lim — - Im 

4-71- r,r'^0 rr' 



QNcMa, 







47r 



lim Im 

--'-+0 



9l{r')gl{r) 



dr" 



(A.9) 



where we assumed ao to be real. 

Let us now analyse the origin of the imaginary part in Eq. (|A.9p . It will be convenient 
to express the r-dependence in terms of the dimensionless variable g = raM, where a = 
g^Cpj^-'K. In these units, the homogeneous Schrodinger equation (Eq. (j5.ip l reads 



^ \ ' +- + 0(1) 



g'Ao) = , 



implying 



9\{q) 



J+i 



1 



-d+' + ... . 



2(/ + l) 

At some order the solution also develops an imaginary part; let us write an ansatz 

1 



(A.IO) 
(A.ll) 



aAo) 



J+2 



+ . . . + i-iiQ , 71 G 



(A.12) 



2(/ + l)" 

The imaginary part in Eq. (j4.2Up behaves as ~ 472 at small q. Inserting into the Schrodinger 
equation, we get for the leading imaginary term 



r — 2r / 



l)-l{l + l)]+il2Q' ■ Q' 







(A.13) 



implying x = / + 5. 

Returning to Eq. ()A.9p , there are in principle two possibilities for the origin of the imaginary 
part. However, according to Eqs. (jXTTj) . (fXT3]) . lim^'^o Imb?/?"'] X.?" dr" Re{l/[ff°(r")]^} ~ 
limr/^o('^')^/(^') — 0- Therefore the imaginary part can only arise from Im{l/[5^(r")]'^}. 
Inserting the asymptotic form of Re[5(J?/r'] from Eq. (jA.lip : using the variable q; and noting 
that this corresponds to the choice oq = aM, we then obtain Eq. (j4.25p . 

It is useful to crosscheck that Eq. ()4.25p produces the correct result in the free limit. In 
the free case there is no i| Iml/>(r)| in Eq. (j4.20p . and a factor ie = iO^ needs to be inserted 
instead, to pick up the correct (retarded) solution. In dimensionless units, the homogeneous 
equation then becomes 

92 Cj' 



+ 



9r{Q) = 



where uj' = uj' /M. We denote 



k = \l ^ + ie . 



(A.14) 
(A.15) 
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The solution with the correct behaviour around the origin (with uq = aM) reads 

g'^{g) = ^sm{kg) . (A. 16) 



We can write 



1 , d f cos[kg') ^ 



The integral in Eq. (j4.25p can now be carried out; the substitution at the upper end gives 
a contribution from the exponentially growing terms exp^—ikg'), present both in the cosine 
and in the sine. Their ratio gives —i, and the total is then 

^-^^ = ^—iui\{ — + ie] i} = -—^9{uj'){Cb')^/^ . (A.18) 



M2 4tt iVa^ y J 27r 

This indeed agrees with Eq. (13.241) . 

A. 2. Scalar channel 

In the scalar channel case, the equations to be solved are (|4.16p . (I4.17p . (j4.20p : the ansatz for 
the solution is in Eq. (I4.22[) . with A given by Eq. (jA.Sp . 

Let us first work out the contribution from the mode / = (S-wave). According to 
Eqs. (liTfll . K22t\ . (fXHl) . the relevant term of denoted by is 

5onu^';r,r') = -^^Ag^ir)g^ir') j dr" . (A.19) 

Inserting into Eq. ()4.16p . making use of Eq. (lA.Sp . and going over into the dimensionless 
variable we get 



M2 47r 



According to Eq. (jA.lip . the first term inside the curly brackets is liuig^Q dg{g^ / g) = —1/2, 
so that we get 



M2 47r 



In principle there are again two possible origins for the imaginary part. However, as we 
saw in the vector channel case, Im[(7[?/f)'] f^dg"Re{l/[g^{g")]'^} ~ (g')'^, so that a non-zero 
contribution can only arise from Im{l/[(7j?(£>")]2}. Furthermore, the derivative can only act 
on the combination multiplying the integral, since 
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Making use of linig/^o d^,' (5!?/^') = —1/2, the S-wave contribution to the scalar spectral 
function thus becomes 



M2 



dg Im 



1 



g0{g)=g-ey2+.. 



(A.23) 



In other words, comparing with Eq. (j4.25p . 6op^{uj') = —a'^p^{uj')/12] the factor is a man- 
ifestation of the suppression ~ V^/M^ apparent in Eq. ()3.3ip . combined with the parametric 
order of magnitude of Vr/M from Eq. (j3.ip . 

Consider then the contribution from the mode I = 1 (P-wave). The relevant term from 
Eq. (|i:22]) . denoted by ^i^-, is 



m=— 1 



Hence we wih need 



— cos 6*, Yi±i( 



Ta/t;— smb'e ^ 
V Svr 



(A.24) 



(A.25) 



In order to take the derivatives in Eq. (j4.16p . we stay with radial coordinates, so that 

^ d 1 d Id , 

Vr = e,— + ee-— +e0— -— . A.26 
or r Of) r sm d 0(p 

Moreover, we choose again r < r', so that r< = r, r> = r' . We will set Q' = 0, after taking 
the derivatives in Eq. (j4.16p . so that the basis is orthogonal. Making use of Eqs. (|A.26p . 
(|A.25p . the terms m = ±1 both yield 



Vr' • Vr 



A- 



d 



Svr [ dr 
while the term m = yields 



_d_ 

dr' 



cos^e + 1 



(A.27) 



Vr' • Vr 



A 



_3_ r 9_ 

47r \ dr 



d_ 

dr' 



giir') 



cos^e + ^^sin^^ 



/\2 



Summing together, we get 
Vr' • Vr ^i^* : 



47r \ dr 



g'Ar) 



d_ 

dr' 



'giirV 



+ 2 



(^0 



gl:{r)gi{r') 

^2 ^^/^2 



^ ^ . (A.28) 



(A.29) 



We now insert c/^(r) = g^{r), gl{r') = -g}{r') f^J'dr" l/[gr{r")? from Eqs. (IXl) . dXH]) . and 
recall that at small r, (r) ^ = {raM)"^ . Thereby 



V w w . ? ,^,3,. r d glig') dg" \ 2gl{g') 

hm Vr'-Vr = - — {aMf hm <^ — / . „,,^ + 7;,V / 



r.r'^O 



dg" 
(A.30) 
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Inserting this into Eq. (j4.16p , and making use of Eq. (|A.8P , we get 



M2 



-a^ lim Im 

27r 



±(9l{Q') 



+ 



[9]{q")? 



(A.31) 



Once again, we need to inspect the origin of the imaginary part. According to Eqs. (jA.lip . 
([XT3|) . Re[gl{Q')] ~ {p'f, Im[c/i (£>')] ~ {p'f, and consequently 



Re 



1 



[9}{q")? 



{e"Y 



Im 



Q'9}{q') 



Im 



{(/r+hiiQ'y 



-HQ 



(A.32) 



so that the only possibility is to consider Im{l/[<^j!^(£)")]2}. The prefactor multiplying this can 
be trivially determined, and we end up with 



M2 



9iVc 
27r 



Q lim / dp Im 



1 

WW 



gi(e)=f2-eV4+.. 



(A.33) 



in analogy with Eq. (j4.25p . Combining Eqs. ()A.23p . (]A.33p . the complete scalar channel 
spectral function can be written as in Eq. (j4.26p . 

To conclude, let us again check that the procedure introduced does yield the correct tree- 
level result. Somewhat unfortunately, the first term in Eq. (j4.26p does not contribute in 
this limit: the subleading term in Eq. (jA.lip would be of 0{q^~^^) in the free case, so that 
9i/q ~ q"^ in Eq. (|A.20p . and 6qp^ vanishes. However, the second term in Eq. (|4.26p survives. 
In dimensionless variables, the homogeneous Schrodinger equation reads 

d^ 



+ 



UJ 



+ ie 



a'- 



(A.34) 



Since there is no imaginary potential, we have had to introduce e = 0^ to pick up the retarded 



solution. The solution normalised to give the desired small- behaviour \g^{^Q) 



9\[q) 



3 



kg 



where k is from Eq. (lA.lSp . We note that 
~s\ii{kg) T -2 



kg 



cos{kg) 



1 d 

k dg 



cos{kg) 



cos{kg) + kgs'm(kg) 
kgcos{kg) — sm(kg) 



g^ + ...] is 
(A.35) 

(A.36) 



whereby 



p%J) 
M2 



- — -a^ lim Im ■ 

27r £i-^oo 



2tt' 



-a Im 



to 



UJ 



+ ie 



cos{kg) + kgsm{kg) 
kgcos{kg) — sm{kg) 



2tt 



(A.37) 



This indeed agrees with Eq. (I3.37p . 
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